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FORTRAN PROGRAM FOR CALCULATING LEADING- AND 


TRAILING-EDGE GEOMETRY OF TURBOMACHINE BLADES 

by Lawrence F. Schumann 

Lewis Research Center and 
U.S . Army Air Mobility R&D Laboratory 

SUMMARY 

A FORTRAN IV program has been written which calculates leading- and 
trailing- edge circle radii, tangency angles on the leading- and trail- 
ing- edge circles, and stagger angle of turbomachinery blade sections 
using only spline points defining the blade surfaces. The program also 
shifts the origin of the blade coordinates to the leading edge of the 
blade* Required input Includes (m, 6) coordinates of a sufficient 
number of spline points to adequately define the two surfaces of the 
blade. Other required input are the radii from the axis of rotation of 
the leading- and trailing- edges. The output from this program may be 
used directly as the geometrical input for a NASA developed program for 
calculating transonic velocities on a b lade- to-b lade stream surface of 
a turbomachine (TSONIC) . The program may be used for axial, radial, and 
mixed flow turbomachine blades. 


INTRODUCTION 


In recent years, several computer programs have been written which 
calculate the velocities on a b lad e-to-b lade surface of a turbomachine 
(ref. 1, 2, 3). These programs require an accurate geometrical description 
of the blade since spline fit curves (ref. 4) are used to define the blade 
surfaces. The inputs required include blade chord, blade stagger, 
leading- and trailing-edge circle radii, angles of tangency on the leading- 
and trailing-edge circles, and intermediate spline points. These inputs 
are usually obtained from measurements made on a physical layout of the 
blade or from a meridional flow analysis program such as MERIDL (ref. 5,6). 
In either case, it is difficult to specify these inputs with the degree 
of accuracy necessary to insure a smooth blade shape without using a 
costly and time consuming trial-and-error procedure with the computer. 

Small errors in any of the input can result in a spline fit through the 
given input spline points in which the first and second derivatives of 
the curve do not vary smoothly. Since the computation of surface velocities 
is highly dependent on the blade surface curvatures which in turn are 
determined from the values of d6/dm and d^e/dm^ obtained from the 
spline fit curve, it is important that the geometry input be specified 
accurately. If the first and second derivatives are not smooth unrealistic 
peaks and valleys in the surface velocity distribution will occur. 
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This report presents a Fortran IV computer program which analytically 
calculates the leading- and trailing-edge circle radii, tangency angles, 
and stagger angle of a turbomachine blade section using an iterative 
procedure. The required input includes enough blade surface (m, 6) coor- 
dinates to adequately define the general shape of the blade surfaces and 
the radius from the axis of rotation of the leading- and t railing-edges . 

This input is much easier to specify than the input for TSONIC. It will 
also yield smoother blade shapes. The output from this program may be 
used directly as the geometrical input for TSONIC (ref. 1) and TANDEM (ref. 3). 
A complete description of the input required and output obtained are given. 

Also given are two example cases and the program listing. 


SYMBOLS 

c meridional chord length, meters, see fig. 1(b) 

m meridional distance, meters, see fig. 1 

r^ radius of leading edge circle, meters, see fig. 1(b) 

r T radius of trailing edge circle, meters, see fig. 1(b) 

R radius from axis of rotation to the blade surface, meters, see fig. 1 

3 angle between the blade surface and the meridional direction, degrees, 

see fig. 1(b) 

6 angular distance between the origin of the input blade coordinates 
and the center of the leading edge circle, radians, see fig. 1(b) 

6 angular coordinate about the axis of rotation, radians, see fig. 1 

0 blade stagger angle from the center of the leading edge circle to 
the center of the trailing edge circle, radians, see fig. 1(b) 

Subscripts : 

L leading edge 

LT leading edge tangency point 

T trailing edge 

TT trailing edge tangency point 

1 surface 1 

2 surface 2 
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DESCRIPTION OF INPUT AND OUTPUT 

The program computes leading- and trailing- edge circle radii, tangency 
anglea, blade stagger and blade chord from input blade surface coordinates. 
The blade surface coordinates are defined by meridional distance, m, and 
the angle about the axis of rotation, 0. Also required as input are the 
radii from the axis of rotation of the leading- and trailing-edges. The 
input blade geometry is shown in figure 1(a). Since spline fit curves are 
used to define the two blade surfaces, the input spline points should not 
reflect the leading- and trailing- edge circle curvatures but should define 
a smooth blade with the surfaces extrapolated to the leading and trailing 
edges as shown in figure 1(a). This is the type of blade definition obtained 
from a known mean camber line and tangential thickness distribution. 

The general procedure for the solution of the problem is described in 
the appendix. An iterative procedure is used which fits a circle in the 
region bounded by the two blade surfaces and the leading or trailing edge 
(fig. 1). The requirement that the circle be tangent to all three boundaries 
leads to a unique solution. Once the radius of this circle is determined, 
the remaining geometry is easily calculated. The variables calculated by 
the program are shown in fignure 1(b). 


INPUT 

The input data sheet along with the values used in an example case 
are shown in figure 2. NP1 and NP2 are integer variables and must be 
right adjusted on 5-column format fields. The remaining input are real 
variables and are specified in 10— column fields. Either English or SI 
units may be used. The input variables are as follows: 

NP1, NP2 Number of spline points on surfaces 1 and 2 respectively 

(maximum 50) . 

Ml, M2 Arrays of m-coordinates of spline points on the two blade 

surfaces, meters, see fig. 1(a). The origin of the m 
coordinates may be anywhere but must be the same for both 
surfaces. 

THETA1 , THETA2 Arrays of 0-coordinates of spline points corresponding 

to Ml and M2, radians, see fig. 1(a). The origin of the 
©-coordinates may be anywhere but must be the same for 
both surfaces. 

ML, MT Meridional distance, m, of the leading and trailing edges 

respectively, meters, see fig. 1(a). The origin for ML 
and MT must be the same as for Ml and M2. 
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RL, RT Radius of blade section from the axis of rotation at the 

leading and trailing edges corresponding to ML and 
MT, meters. 

Surface 1 must always be the upper surface, i.e., the surface with 
the larger n ume rical values of 0 as shown in figure 1(a). If this is 
not the case, an error message will be printed and no further computation 
will be done. Also, the leading- and trailing edge thicknesses must not 
equal zero for this technique. 

It is advisable, if possible, to specify spline points beyond the 
leading and trailing edges of the blade as shown in figure 1(a) since this 
reduces the effect of the end conditions used in the splint subroutine 
(i.e. that the second derivative at the end points is one-half of that at 
the adjacent point). However, this is not an absolute requirement since 
the splint subroutine has the capability of linear extrapolation and there 
fore the specified spline points do not have to extend up to or beyond 
the leading- and trailing-edges. 


OUTPUT 

Two examples of the output obtained from the program are shown in 
figure 3. The output is divided into three points. The first part is a 

listing of the input data. The second part is a listing of the variables 

computed by the program. They are: 

RLE r^, meters, see figure 1(b) 

BETAL1 ^lxi* degrees , see figure 1(b) 

BETAL2 &lt 2* degrees » see fi 8 ure 

RTE r T> meters, see figure 1(b) 

BETAT1 e xTl' degrees , see figure 1(b) 

BET AT 2 $ tt2 > degrees, see figure 1(b) 

CHORD c, meters, see figure 1(b) 

STGR a, radians, see figure 1(b) 

DELTA 6, radians, see figure 1(b) 

The third part of the output is a list of the input blade coordinates 

with the origin shifted to the leading edge of the blade as shown in figure 1(b). 

Also shown are the calculated tangency points. The values of the first and 
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second derivatives (d0/dm, d 6/dm ) obtained from the spline fit are 
also listed. A check on the smoothness of the blade can be made by 
examining the values of the second derivatives. They should vary smoothly 
from inlet to exit and have no wild oscillations. Any input coordinates 
whose m value is less than m^ or greater than mpx are not listed 
since these values should not be used as input to TSONIC. The values of 
the tangency points also should not be used as input coordinates to TSONIC 
(except as dummy variables) since these values will be calculated by 
TSONIC. The remaining coordinates ma> be used directly as input to 
TSONIC with no further modification. 


APPLICATION OF THE PROGRAM 

The primary purpose of this program is to easily obtain the input for 
TSONIC or TANDEM when only the blade surface coordinates are known or can 
easily be obtained (e.g. from a known mean camber line and tangential 
thickness distribution - fig. 4(a)). This is the type of output obtained 
from a meridional flow analysis program such as MERIDL. The current version 
of MERIDL (ref. 5 and 6) employs the technique presented in this report 
for the calculation of TSONIC input data. This program can also be used 
in conjunction with a physical layout since it is often difficult to 
measure angles and leading and trailing edge radii with any degree of 
accuracy while it is a relatively easy task to obtain 5 or 6 spline points 
on the blade surfaces which will yield a smooth blade shape (fig. 4(b))* 

It is also easy to obtain splitter blade geometry using this program by 
simply changing the value of m^ and Rl as is shown in example 2 
(figs. 3(b) and 4(c)). 



PROGRAM LISTING 


OOOC 100 
QC00200 
OC 00 300 
0000400 C 
0C0050C C 
OCOC600 C 
0000700 
0000900 
0000900 
0C01000 
0001100 
0001200 
0001300 
0001400 
- 001 500 
000160C 
00017*0 
OC 0 1 80 r 
0 0 0 1 9 r 0 
0002000 
0002100 
0002200 
CO 02 300 
C00240C 
0C025C0 
CC02600 
0002700 
0002800 C 
0002900 C 
C0O3000 C 
0003100 
OC 03200 
0003300 
0CC3400 
0C03500 
OC 03600 
OC 03700 
OC 03800 
0003900 
OCX 4001 
0004100 
0004200 
0004300 
0004400 
0004500 
CC046C0 
0004700 
0004 800 
0004900 
0005000 
0005100 
0005200 
000530C 
0005400 


DIMENSION Ml (50) ,M2 (5C) , THETA 1 (50) , THFTA2 (50) , DTDM1 (50) 
1D2TDM1 (50) , D2TDH2 (50) 

REAL Ml, M2, ML, NT,MTL1,MTL2, MTT1,MTT2 

PFAD AND WRITE INPOT DATA 

10 PEAC(5, 1000, END*999) 

WRITE (6,1001) 

WRITE (6, 1000) 

WRITE (6,1100) 

PEAC(5, 1010) NP1, NP2 
WRITE (6,1010) NP1 ,NP2 
WRITE(6, 1110) 

READ (5, 1020) (Ml (I) ,1*1, NP1) 

WRITE (6, 10 30) (Ml (I) ,I*1,NP1) 

WRITE (6,1120) 

READ (5, 1020) (THETA 1 (I) ,1*1, NP1) 

WRITE (6,1030) (THETA 1 (I) ,I=1,NF1) 

WRITE(6, 1130) 

READ (5 , 1 020) (M2 (I) , I* 1 , N P2) 

WRITE(6, 1030) (M2 (I) ,I*1,NP2) 

WRITE (6,1140) 

READ (5, 1020) (THETA 2 (T) ,1*1, NP2) 

WRITE (6,10 30) (THETA 2 (I) , I s 1 , NP2) 

WRITE(6, 1160) 

READ (5 , 1 020) ML, MT, R L, RT 
WRITE (6, 1030) ML,HT,RL,RT 

COMPUTE LEADING EDGE RADIUS AND TANGENCY ANGLES 

CALL SPLINT (Ml , THETA 1 ,N PI, ML, 1,THL1, DL1 ,D2L1) 

CALL SPLENT ( Ml , NEl ,THET A1 , DTDMl , D2TDM1) 

CALL SPLINT (M2, THETA2, NP2, ML, 1,THL2, DL2,D2L2) 

CALL SPLENT ( M2, NP2, THETA 2, DTDH2, D2TDM2) 

DA HP* 1 • 

ITER=0 

20 BETAL1*ATAN (RL*DL1) 

BET AL2* ATAN (RL*DL2) 

RLE 1* (THL1-THL2) *RL*COS( (BETAL1 ♦BETAL2) /2 , ) /2 . 

IF (RLE1 .LE.O.) GO TO 120 
30 HTL1*HL»RLE1*(1.-SIN(BETAL1)) 

MTL2 = HL» BLE1 * (1 . *SI N (BET AL2) ) 

CALL SPLINT (Ml, THETA 1,N PI, HTL1, 1,TTL1,DTL1 ,D2TL1) 

CALL SPLINT (M2 ,THETA2,NP2, MTL2, 1,TTL2, DTL2, D2TL2) 

BET AL 1 * ATAN (RL* DTL1 ) 

BETAL2*ATAN (RL*DTL2) 

RLE2*RL* (TTL1-TTL2)/ (COS (BETALl) *005 (BET AL2) ) 
IP(ABS((RLE2-RLE1)/RLFl) .LT..C001)GO TO 50 
ITEB*ITER* 1 

IP (ITER. LE. 100) GO TO 40 
WRITE (6,1200) 

GO TO 50 

40 RLE1*(DAHP*RLE1+FLE2)/(DAMP+1») 

IP ( RL E 1 • GT • 0 . ) GO TO 30 


DTDM2 ( c *) 
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00055C0 
CG0560C 
0005700 C 
000580C C 
0005900 C 
0C06000 
0006100 
0006200 
0006300 
0006400 
000650C 
0006600 
0006700 
0006800 
0006900 
0007000 
0007100 
0007200 
0007300 
0007400 
0007500 
0007600 
0C07700 
0007800 
0007900 
0008000 
0008100 
0008200 
0008300 
0008400 C 
0008500 C 
0008600 C 
0008700 
0008800 
0008900 
0009000 
0009100 
0009200 
0009300 
0009400 
0009500 
0009600 
0009700 
0009800 
0009900 
0010000 
0010100 
0010200 
0010300 
0010400 
0010500 
0010600 
0010700 
0010800 


DAflP*DAHP+ 20 • 

GO TO 20 

COMPUTE TRAILING EDGE RADIUS AND TANGENCY ANGLES 

5C CALL SPLINT (Ml # THETA 1 ,NPl,RT,1,THT1,DT1,D2T1) 

CALL SPLINT (M2, THET A2 ,NP2, MT, 1 , THT2 , DT2,D2T2) 

DAMP* 1 • 

ITER*0 

60 DETAT1 * ATAN (RT*DT1 ) 

BETAT2*ATAN (RT*DT2) 

RTE1* (THT1-THT2) *RT*COS ( (BETATl +BETAT2) /2. ) /2. 

IF(RTE1.LE.O.)GO TO 120 
70 MTT1*MT-RTE1* (1.+SIN (BETATl) ) 

MTT2*HT-RTE1* (1 . -SIN (BETAT2) ) 

CALL SPLINT (Ml, THET A 1 ,NPl , MTT1 , 1 ,TTT 1, DTT 1, C2TT i) 

CALL SPLINT (M2, THET A2,NP2,HTT2, 1 t ,TTT2,DT72, D2TT2) 

BETATl * ATAN (RT*DTT1 ) 

BET AT2*ATAN (RT*DTT2) 

RTE2*RT* (TTT1-TTT2) / (COS (BETATl) ♦COS (BETAT2) ) 
IP(ABS((RTE2-RTEl)/RTEl).LT..0001)GO TO 90 
ITER* ITER* 1 

IF (ITER.LE.100) GC TO 80 
WRITE (6, 1210) 

GO TO 90 

80 RTE1*(DAMP*RTE1*RTE2)/(DAHP*1.) 

IF(RTE1.GT.0.)GO TO 70 
DAHP*DAHP*20 • 

GO TO 60 

CALCULATE AND PRIUT OUTPUT DATA 

90 DELTA* (TTL1*C0S (BET AL2) ♦TTL2*COS (BETAL1) ) / (COS ( BET AL 1) ♦COS (BETAL2) 

D 

STGR* (TTT1 *COS (BETAT2) ♦TTT2*COS (BET ATI) )/ (COS (BETATl) ♦COS (BETAT2) ) 
1-DELTA 
CHORD*HT-HL 
BETAL1*BETAL 1*57. 29578 
BETAL2*BETAL2*57. 29578 
BET ATI *BITAT 1*57. 29578 
BETAT2*BETAT2*57. 29578 
MTL 1*MTL 1-ML 
MTL2* HTL2-ML 
MTT1*HTT1-ML 
HTT2*HTT2- ML 
TTL1*TTL1-DELT A 
TTL2*TTL2-DELTA 
TTT1*TTT1- DELTA 
TTT2*TTT2- DELTA 
WRITE(6, 1160) 

WRITE (6,1030) RLE1,BETAL1 , BET AL2, RTE1 , BET ATI , BETAT2 
WRITE (6, 1170) 

WRITE (6, 1030)CHORD, STGR, DELTA 
WRITE (6,1180) 
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0010900 
0011000 
0011100 
0011200 
0011300 
0011400 
0011500 
C0 11600 
0011700 
0011800 
0011900 
0012000 
0012100 
OC 12200 
0012300 
00 1240 0 
00125CC 
0012600 
0012700 
00128C0 
0012900 
001300^ 
00 13 IOC 
0013200 
0013300 
0013400 
0013500 
0013600 
0013700 
0013800 
0013900 
001400C 
0014100 
0014200 
0014300 
0014400 
OC 14500 
0014600 
0014700 
0014800 
0014900 
0015000 
0015100 
0015200 
0015300 
0015400 
0015500 
0015600 
0015700 
0015800 
0015900 
0016000 
0016100 
0016200 


J«1 

WRITE (6 f 1040) J, HTL1 ,TTL1 ,DT Li , D2TL1 

DO 100 1*1, NP1 

Hi (I) *Hl (I) -HI 

THETA 1(1) *THETA 1 (I) -DELTA 

IF (Hi (I) .LE. HTL1.0R.Hl (I) .GE. HTT1)G0 TO 100 
J*JO 

WRITE (6, 1050) J, HI (I) ,THETAl (I),DTDH1(I) , D2TDH 1 (I) 

100 CONTIWO E 
J*J4l 

WRITE (6, 1040) J, HTT1 ,TTT1 ,DTT1 ,D2TT1 
WRITE (6, 1190) 

J«1 

WRITE (6 , 1040) J, HTL2,TTL2,DTL2, D2TL2 

DO 110 1*1, HP2 

H2(I)*H2(I)-HL 

THETA2 (I) *THETA2 (I) -DELTA 

IF (H2 (I) . LE. HTL2.0R, H2 (I) . GE. HTT2) GO TO 110 
J*J*1 

WRITE (6,1050) J, H2 (I) ,THETA2 (I) ,DTDH2(I) ,D2TDH2(I) 

110 CONTINUE 
J* J ♦ 1 

WRITE (6, 1040) J, HTT2,TTT2,DTT2,D2TT2 
GO TO 10 

120 WRITE (6, 1220) 

GO TO 10 

i'ORH AT STATEHENTS 

1000 FORHAT (80H 

1 ) 

1001 FORHAT (1H1 ) 

1010 FORHAT (1615) 

1020 FORHAT (8F10. 6) 

1030 FORHAT (8G1 5. 6) 

1040 FORHAT (6X, 1 4HTANGENCY POINT , 3X, 12 , 4G 15. 6) 

105C FORHAT (23X,I2,4G15.6) 

1100 FORHAT (/10H NP1 NP2) 

1110 FORHAT (8X, 8HH1 ARRAY) 

1120 FORHAT (81, 12HTH3TA1 ARRAY) 

1130 FORHAT (8X, 8HH2 ARRAY) 

1140 FORHAT (8 X, 1 2HTHETA2 ARRAY) 

1150 FORHAT (6X,2HHL, 13X,2HHT, 13X ,2HRl, 13X, 2HRT) 

1160 FORHAT (///32X,25H** COHPOTED ELADE DATA ••///6I, 3HRLE, 1 21, 
16HBETAL1,9X, 6HBETAL2, 9X , 3HRTE, 1 2X,6HPETAT1 , 9X ,6HBETAT2) 

1170 FORHAT (6 X,5HCHORD, 1 0X ,4HSTG R, 1 1 X, 5H DELTA) 

1180 FORHAT (///34X,2 1 H* BLADE COORDINATES */30X, 30 H (ORT GIN AT BLADE LEA 
1DING ED' E) //40X , 9HSUR FACE 1//24X, 1HI,6X, 1 HR, 14X,5H THETA, 10X,4HDTDH 
1 , 1 1 X, 6HD2TDH2) 

1190 FORHAT (//40X , 9HSURFACE 2//24X ,1HI,6X,1HH, 14X,5HTHETA, 10 X,4HDTDH, 
111 X,6HD2TDH2) 

1200 FORHAT (//5X,61HS0LUTI0N FOR LEADING EDGE HAS NOT CONfERGED IN 100 
1 ITERATIONS) 

1210 FORHAT (//5X,62HSOLUTION FOR TRAILING EDGE HAS NOT CONfERGED IN 100 
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0016300 
0016400 
0016500 
0016600 
0016700 
0016600 
0C16900 
CC 17000 


1 ITERATIONS) 

1220 FORHAT (///10 X # 43HINPUT ERROR — ONE OF THREE POSSIBLE CAUSES/ 

1 17X # 32H1 • SURFACES 1 AND 2 ARE INVERTED/ 

2 1 7X#54H2 • BLADE THICKNESS IS ZERO AT LEADING OR TRAILING EDGE/ 
317X # 68H3, SURFACES 1 AND 2 CROSS OVER AT SORE POINT DOE TO A BAD S 
4PLINE FIT) 

999 STOP 
END 


onoo ioo 

0000200 C 
C000300 C 
0000400 C 
0000500 C 
C000600 C 
0000700 C 
0C00900 
0001000 
0001100 
0001200 
0001300 
0001400 
0001500 
0C01600 
0001700 
0001800 
0001900 
0002000 
0C021C0 
0002200 
0002300 
0002400 
0002500 
0002600 
0002700 
0002800 
0002900 
0C03000 
0003100 
0003200 
0003300 
0003400 
0003500 
C003600 C 
0003700 
0003800 
0003900 
0004000 
0004100 
0004200 
0004300 
0004400 
0004500 
0004600 
0004700 
0004800 
C.004900 
0005000 
0005100 
0005200 
0005300 
0005400 
0005500 


SUBROUTINE SPLINT ( X, Y, N , Z, RAX , TI NT , DTD X, 02 YDX2) 

--SPLINT CALCULATES INTERPOLATED POINTS AND DERIVATIVES 
--POR A SPLINE CURVE 

--END CONDITION - SECOND DERIVATIVES AT END POINTS APE 
--SDR1 AND SDRN TINES SECOND DERIVATIVES AT ADJACENT POINTS 

DIMENSION X (N) ,Y(N) ,Z (MAX) , TINT (MAX) f DYDX (MAX) , D2TDX2(MAX) 
DIMENSION G (101) , SB (101) , EH (101) 

IERR * 0 
SDR 1 * .5 
SDRN * .5 

TOLER* ABS (X (N) -X ( 1 )) /FLOAT (N) * 1 . E-5 

C * X(2)-X(M 

IF (C. EQ.O.) GO TO 110 

SB (1) * -SDR 1 

G ( 1 ) * 0. 

NO * N-1 

IF (NO.LE.0) GO TO 140 
IF (NO.EQ.1) GO TO 20 
DO 10 1*2, NO 
A * C 

C * X(I*1)-X(X) 

IF (A*C.EQ.O.) GO TO 130 
IF (A*C.LT.0.) IERR * 1 
W * 2 • * (A^C) -A*SB (I -1 ) 

SB (I) * C/B 

F * (V (IM)-Y(I) )/C-(Y(I) - Y (1-1 ) ) /A 
10 G ( I ) * (6.*F-A*G (1-1) )/B 
20 EH ( N) * SDRN*G (N-1) /(I. ♦SDRN*SB (N-1) ) 

DO 30 1*2, N 
K * N ♦ 1 - 1 

30 EH (K) * G(K) -SB(K) *EH(KO) 

IF (HAX.LE.O) RETURN 

ENTRY SPLENT (Z , HA X , Y INT , DYDX , D2 YDX2) 

DO 120 1*1, MAX 
K*2 

IF (ABS (Z (I) -X (1)) . LT.TOLER) GO TO 40 
IF (Z(I) ,GT.2.0»X(1j-X(2)) GO TO 50 
GO TO 80 

40 TINT (I) * Y ( 1 ) 

SK * X (K) — X (K— 1) 

GO TO 110 

50 IF (ABS(Z(I)-X(K)) .LT.TOLER) GO TO 60 
IF <Z(I) .GT.X(K) ) GO TO 70 
GO TO 100 
60 YINT(I) * Y ( X) 

SK * X (K) -X (K-1) 

GO TO 110 

70 IF (K.GE.N) GO TO 90 
K * 

GO TO 50 

80 S2 * X (2) -X ( 1 ) 
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3005600 
0005700 
0005800 
0C0590C 
0C 06000 
0006100 
0006200 
0006300 
0006400 
00065C0 
0006600 
0006700 
0006800 
0006900 
0007000 
0007100 
OC 07200 
0007300 
0007400 
CC07500 
0007600 
0007700 
00078GC 
000790f 
0008000 
0008100 
0008200 
0008300 
0008400 
0C085C0 
0008600 
0008700 
0009800 
OC 08900 
0009000 
0009100 
0009200 
0009300 
0009400 


YO * EH(1) •S2**2*2.*Y(1)-Y(2) 

DYDX(I) * (Y(2)-Y(1))/S2-7.*EH(1)/6.*S2 
YINT(I) * Y0*DYDX(I) *(Z(I) -X ( 1 1 ♦S2) 

D2YDX2 (I) » 0. 

00 TO 120 

90 IP (Z(I) ,LT.2.*X(N) -X (N-1)) GOTO 100 
SN ■ X <N) -X (W- 1 ) 

YNP 1 * EPI(N) *SN**2*2. *Y <N) -Y (N-1) 

OYDX(I) * (Y (N) -Y(N-1))/SN*7.*EH(N)/6.*SM 
YINT(I) = YNP1*DYDX(I)*(Z(T) -X (N) -SN) 

D2 YDX2 (I) * 0. 

GO TO 120 

100 SK * X(K)-X(K-1) 

YINT(I) * EH(K-1)*(X(K) -Z(I) ) **3/6./SK *EH (K) * .2 (I) -X (K-1 ) ) •• V6. 

1 /SK+(Yr<)/SK -EN(K)*SR /6.) •(Z(I)*X(K-1)) *(Y(K-1)/SK -FN(K-I) 

2 * SK/6 • # * (X (K) -Z(I) ) 

110 DYDX(I) »-E«(K-1) *(X (K) -Z (T) ) **2/2.0/SIC ♦ FB <K> M* (K-1) -Z (I) ) **2/2. 

1 /SK* (Y(K) -Y (K-1 ) )/SK - (Efl (K) -EH (K-1) ) *SK/6. 

D2 Y0X2 ( I) * EH(K)-(X(K)-Z(I))/SKMEH(K)-En(K-1)) 

120 CONTINUE 

IP (IERR.EQ.O) RETURN 
UO WRITE (6, 1000) 

WRITE (6 , 1020) N, (X(I) ,Y(I) ,I*1,N) 

IP (IERR.EQ.O) STOP 
WRITE (6 , 1030) 

RETURN 

140 WRITE(6, 1010) 

WRITE (6 , 1020) N, (X(I) ,Y (I) ,1-1, N) 

STOP 

1000 FORHAT { 1H 1 , 10 X , 44 HS PLI NT ERROR — ONE OP THREE POSSIBLE CAUSES/ 

1 1 7X , 5 1 HI • ADJACENT X POINTS ARE DUPLICATES OP EACH OTHER./ 
217X,38H2. SORE X POINTS ARE OUT OP SEQUENCE. / 

3 1 7x , 32 H3 • SORE X POINTS ARE UNDEPINED.) 

1010 FORHAT ( 1 H 1 , 10 X , 62 HS PL I NT ERROR — NHHBER OF SPLINE POINTS GIVEN I 
IS LESS THAN TWO) 

1020 FORHAT (//17 X, 1 fl HNU HBER OP POINTS * , I4//1 7X,8 HX ARRAY, 6l # 8HY APR 
1AY/(17X,2G13.5) ) 

1030 FORHAT (1H1) 

END 
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APPENDIX 

CALCULATION PROCEDURE OF THE PROGRAM 

The program computes the leading- and t railing-edge circle radii and 
tangency angles by an iterative procedure. Figure 1 shows the basic blade 
geometry. The input coordinates are given as (m,6) coordinates where m 
is the length in the meridional direction and 6 is the angle in radians 
about the axis of rotation. The origin for both the m and 6 coordinates 
may be anywhere since it will be shifted by the program to the leading edge 
of the blade. 

From figure 1(b) (with the assumption that rT << R^) , the following 
equations can be written for the leading edge circle: 


°LT1 " m L + r L (1 ‘ 8la 

S LT1 ) 

(1) 

“LT2 " “L + r L (1 + 8i “ 

B LT2^ 

(2) 

_ r L c °s S m 

9 lti \ 


(3) 

. r L C0s S LT2 
®LT2 ' 6 ‘ ^ 


(4) 


Equations (3) and (4) can be solved for r which gives 


_ *L (9 LT1 ~ 9 LT2 ) 
r L ’ cos e LT1 + cos e LT2 

The tangency angles are determined from the equations 
tan 6 LT1 

tan B LT2 “ ®L Gta ) LT2 


(7) 
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The angle 6 and d0/dm for a given value of m are obtained from a 
spline fit curve (subroutine SPLINT, ref. 4) on the input surface coordinates. 
Equations (1), (2), (5), (6) and (7) along with the spline fit curves. on 
the surfaces are sufficient to solve the problem. Similar equations can 
be written for the trailing edge. 


The program procedure is as follows: 


(1) A first estimate of the leading edge radius is obtained from the equation 


r L a* *2 (6 U ‘ ®L2 ) co8 



( 8 ) 


where 6 and 0. are values obtained from a spline fit curve at the leading 
edge (see fig. 1(a)). 

(2) Values of and m LT 2 are obtained from equation (1) and (2) 

using the e - 'mated value of r^, the leading edge values of and the 

input value rY 

(3) Splint calls are done on both surfaces using the computed value.: ^ 

and m LT 2 which will yield values of 0^ and B^* 

(4) A new value of leading edge radius is obtained from equation (5). 

(5) An average value of r L is used to start the next iteration. The average 
value will be weighted toward the previous value if r-^ becomes negative. 

The values of &LT from the previous iteration are used. 

(6) When the values of r^ for two successive iterations are within 0.01% 

of each other, the solution is converged. 

(7) The solution for the trailing edge is obtained. 

(8) Values of 6, stagger (a), and chord (c) are obtained from the equations 


9 LT1 008 6 LT2 * 9 LT2 COa g LTl 
cos 8 m + cos S LT2 


(9) 


0 mm - cos fi. iU lir , + cos 
TT1 TT2 TT2 TT1 

C0S ^TTl + C0S 6 TT2 


- 6 


( 10 ) 
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0 

1 

hF 

1 

(ID 

(9) The Input blade coordinates are shifted so that 
blade leading edge by the following transformation 

the origin is at 

m - m - 

(12) 

CD 

■ 

CD 

1 

o> 

(13) 
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Figure 3(a) EXAMPLE OUTPUT 

EXAMPLE 1 - CENTRI PUGAL COMPRESSOR HUB SECTION 

NP 1 NP2 
6 6 

HI ARRAY 


0.000000 

0 • 257480 E-0 1 

0 • 852690E-0 1 

C. 149890 

*.207820 

*. 23172° 

THETA 1 

ARRAY 





0.504740 

0.833650 

1.10168 

1.17271 

1.25245 

1 . 29961 

M2 ARRAY 





o.ooococ 

0. 257480F-01 

C.852690E-01 

0. 149890 

0.2*7820 

*.23172° 

THETA2 

ARRAY 





0.420240 

0.76123C 

1.05180 

1.13873 

1.22981 

1.2813* 

ML 

MT 

RL 

RT 



C. 000000 

0.231720 

0 , 833300E-01 

0.264500 




COMPUTED BLADE DATA •• 


RLE 

BETAL1 

3FTAL2 

RTF 

BET ATI 

0.217386E-02 

51.2904 

51.2865 

0.213013E-02 

28. 8679 

CHORD 

STGR 

DELTA 



0.231720 

0.790331 

0.495592 




BF^AT 
3 1. 1041 


♦ BLADE COORDINATES • 
(ORIGIN AT BLADE LEADING EDGE) 

SURPACE 1 




I 

M 

THETA 

DTD M 

TANGENCY 

POINT 

1 

0 • 477540E-03 

0.1631 35E-C1 

14.9739 



2 

0 • 257480E-0 1 

0.338058 

9.94511 



3 

0 . 852690E-01 

0.606088 

1.46705 



4 

0.14989* 

0.677117 

1.04757 



5 

0.207820 

0.756858 

1.79708 

TANGENCY 

POINT 

6 

0.228561 

0.797385 

2.08430 


D2TDM? 

-134.296 
-263.7*2 
-21. 1741 
8. 19146 
17. 6849 
10 . 011 ? 


SURFACE 2 




I 

M 

THETA 

DTD M 

TANGENCY 

POINT 

1 

0.387008E-02 

-0,1631 5 3E-01 

14.9718 



2 

0 • 257480E-01 

0.265638 

10.3914 



3 

0 • 852690E-01 

0.5562C8 

1.77377 



4 

0. 149890 

0.643137 

1.25508 



5 

0.207820 

0.734218 

1.98404 

TANGENCY 

POINT 

6 

0.230693 

0.783442 

2. 28825 


D2TDM2 
-152. R93 
-266.830 
-23. 7373 
7.68391 
17.4830 
9.11722 
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Figure 3(b). - EXAMPLE OUTPUT 

EXAMPLE 2 - CENTRIFUGAL COMPRESSOR HUB SECTION - SPLITTER BLADF 

NP 1 NP2 
6 6 

Ml ARRAY 


0.000000 0 • 257480E-0 1 

THETA1 ARRAY 

0.852690E-01 

0.149890 

0.207820 

0. 231720 

0.504740 0.833650 

M2 ARRAY 

1.10168 

1.17271 

1.25245 

1.29961 

0.000000 0.25748CE-01 

THETA2 ARRAY 

0. 852690E-01 

0.149890 

0.207820 

0. 231720 

0.420240 0.761230 

ML NT 

0.5173 20 E-0 1 0.231720 

1.05180 

RL 

0.102740 

1.13873 

RT 

0.264500 

1.22981 

1.2B139 




*• COMPUTED 

BLADE DATA •• 



RLE 

BETAL1 

BETAL2 

RTE 

BETAT1 

BFTAT2 

0.284874E-02 

23. 3568 

23. 5919 

0,21301 3E-02 

28.8679 

31.1841 

CHORD 

STGR 

DELTA 




0.179988 

0.288974 

0.996949 





* BLADE COORDINATES • 
(ORIGTN AT BLADE LEADING EDGE) 

SURFACE 1 




I 

N 

THETA 

DTD M 

D2TDM2 

TANGENCY 

POINT 

1 

0.171 933E-02 

0.254581 E- 01 

4.20327 

-150.820 



2 

0. 335370E-01 

0.104730 

1. 46705 

-21. 1741 



3 

0.98 1 579E- 01 

0.175760 

1.04757 

8.1<*146 



4 

0.156088 

0.255501 

1.79708 

17. 6848 

TANGENCY 

POINT 

5 

0.176829 

0.296028 

2.08430 

10.M1O 






SURFACE 2 





I 

H 

THETA 

DTDH 

D2TDH2 

TANGENCY 

POINT 

1 

0 • 398 887E-02 

-0 . 2541 26F-0 1 

4. 25075 

-143.920 



2 

0. 33537 OE -01 

0 . 548503E-01 

1.77377 

-23.7373 



3 

0.98 1 579E-01 

0.141780 

1.25508 

7.68391 



4 

0.156088 

0.232860 

1.98404 

17.4830 

TANGENCY 

POINT 

5 

0. 178961 

0.282084 

2.28825 

9.11722 
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